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We determine the commutation relations satisfied by the quantized electromagnetic field in the 
presence of macroscopic dielectrics and conductors, with arbitrary dispersive and dissipative prop- 
erties. We consider in detail the case of two plane-parallel material slabs, separated by an empty 
gap, and we show that at all points in the empty region between the slabs, including their surfaces, 
the electromagnetic fields always satisfy free-field canonical equal-time commutation relations. This 
result is a consequence of general analyticity and fall-off properties at large frequencies satisfied by 
the reflection coefficients of all real materials. It is also shown that this result does not obtain in the 
case of conductors, if the latter are modelled as perfect mirrors. In such a case, the free-field form of 
the commutation relations is recovered only at large distances from the mirrors, in agreement with 
the findings of previous authors. Failure of perfect-mirror boundary conditions to reproduce the 
correct form of the commutation relations near the surfaces of the conductors, suggests that caution 
should be used when these idealized boundary conditions are used in investigations of proximity 
phenomena originating from the quantized electromagnetic field, like the Casimir effect. 
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I. INTRODUCTION 



The interaction of radiation with matter has always 
been a fascinating subject of investigation, and in fact it 
is at the roots of quantum mechanics, with Planck's work 
on black body radiation. Even though, after the develop- 
ment of Quantum Electrodynamics (QED) in the middle 
years of last century, all fundamental principles involved 
in this interaction are undoubtedly well understood at 
the microscopic level, recent experimental advances have 
prompted much interest in theoretical studies of the 
quantized electromagnetic (e.m.) field in close proxim- 
ity to macroscopic bodies. A thorough understanding 
of this problem is indeed needed for a correct interpreta- 
tion of numerous important proximity phenomena of e.m. 
origin, that include cavity QED the Casimir effect 
0, radiative heat transfer quantum friction the 
Casimir-Polder interaction of Bose-Einstein condensates 
with a substrate Q, etc. Apart from the intrinsic inter- 
est of these phenomena, it has been shown recently that 
the quantum fluctuations of the e.m. field surrounding 
macroscopic bodies, that are at the origin of the Casimir 
effect, could have exciting application in nanotechnology 
i- 

The common feature of the above e.m. phenomena, 
is that they all involve several macroscopic bodies and 
possibly one or more microscopic objects (atoms, ions 
etc.) placed in a vacuum and separated by distances 
(typical separations range from a few tens of nanometers 
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to several microns) that, while small from a macroscopic 
point of view, are still large compared to the interatomic 
distance in condensed bodies. In such circumstances, the 
microscopic point of view is not of great help, because the 
long range character of the e.m. field implies that macro- 
scopically large number of atoms are inevitably involved 
in the interaction. A much more effective approach would 
be to describe the influence of the macroscopic bodies 
on the quantized e.m. field in the vacuum just outside 
their boundaries, in terms of macroscopic features of the 
bodies like the electric and/or magnetic permittivities. 
On physical grounds, one expects that such an approach 
should be feasible, in certain circumstances at least, be- 
cause the wavelengths of the e.m. fields participating 
in these phenomena are expected to be of the order of 
the bodies separations, and are therefore large on the 
atomic scale. This being the case, use of macroscopic re- 
sponse functions of the bodies should be legitimate. An 
inevitable complication that one faces though, when deal- 
ing with macroscopic response functions of real bodies, is 
that they always display dispersion and absorption. As 
is it well known, the former feature is mathematically 
reflected in the fact that response functions depend on 
the frequency u (we shall neglect spatial dispersion, and 
therefore we shall not consider the possible dependence 
of the response functions on the wave- vector k), while 
the presence of dissipation entails that the response func- 
tions have a non- vanishing imaginary part. The existence 
of absorption, in particular, greatly complicates explicit 
quantization of the macroscopic e.m. field. Unfortu- 
nately, such a difficulty cannot be disposed of by simply 
neglecting dissipation, because dispersive, real- valued re- 
sponse functions inevitably violate causality, and must 
therefore be rejected. 
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Fortunately, though, there exists a way out that avoids 
the above mentioned difhculties. This is so because a full 
quantization of the e.m. field is usually not needed, as 
the quantities of interest are typically statistical averages 
of quadratic expressions involving the macroscopic e.m. 
field. For systems that are in thermodynamic equilib- 
rium, such averages can be expressed in terms of (the 
imaginary part of) suitable macroscopic response func- 
tions, as a result of general fluctuation-dissipation theo- 
rems derived in the framework of linear-response theory 
Q. This general approach was probably pioneered by 
Rytov Q in his investigations of e.m. fluctuations in the 
presence of macroscopic bodies in thermal equilibrium, 
and it was later used by Lifshitz [§| in his famous theory 
of dispersion forces between macroscopic condensed bod- 
ies. In one form or another, the fluctuation-dissipation 
theorem is used in all existing approaches to problems 
involving the quantized e.m. field in the vicinity of or in- 
side macroscopic bodies. In the seventies of last century, 
Agarwal used it as the basis of a systematic investiga- 
tion of QED in the presence of dielectrics and conductors 
[To| . For a review of the most recent work we address 
the reader to Refs.pll. (see also Refs. therein). It 
is important to note that this approach is not restricted 
to systems in global thermodynamic equilibrium, as it 
is still valid in systems that are only in local thermo- 
dynamic equilibrium. This feature permits to include 
within the scope of the theory other important phenom- 
ena, like radiative heat transfer between closely separated 
bodies (for a recent review see @), and quantum friction 
Recently, the theory has also been applied to the 
investigation of Casimir-Polder 0, and Casimir 
forces out of thermal equilibrium. 

In this paper, we reexamine the basic quantum-field- 
theoretical problem of the commutation relations satis- 
fied by the quantized e.m. field in the presence of di- 
electrics and/or conductors, in the framework of the gen- 
eral macroscopic theory described above. Our interest 
in this problem arose from a paper by Milonni [lH on 
the Casimir effect, in which it was found that near a 
perfectly reflecting slab, the transverse vector potential 
and the electric field satisfy a set of equal-time canoni- 
cal commutation relations of a different form form those 
holding for free fields. This result is quite worrisome, in 
view of the very fundamental character of commutation 
relations, because it contradicts one's expectations based 
on microscopic theory, and therefore it deserves detailed 
investigation. We remark that unexpected commutation 
relations between the annihilation and creation operators 
for the e.m. field inside a cavity were also found more 
recently in Ref. [l^. A par tial resolution of the para- 
dox was offered in Ref. [17| . in which the problem of 
the e.m. commutation relations was investigated within 
a simplified form of QED, in one space dimension. By re- 
lying on a simple quantum theory of the one-dimensional 
lossy beam splitter, along the lines of Ref. p"8|_fsee also 
[l9|). it was shown that the anomalies found in 16] in the 
commutators of the annihilation and creation operators 



were associated with a particular choice of the the cavity 
e.m. modes. The authors of Refs. [13] also showed that 
the canonical one-dimensional commutation relations in- 
volving the vector potential and the electric field do not 
display any anomalous behavior. No detailed explanation 
was however provided for the modified form of the equal- 
time commutators derived in Ref. [T5| . apart from the 
remark that the boundary conditions (b.c.) satisfied by 
the e.m. field in the case of ideal mirrors are incompatible 
with the transverse delta function form of the full canon- 
ical commutator in three dimensions. These authors fur- 
ther conjecture that standard equal-time commutators 
would probably be restored after incorporating, in the 
full three-dimensional setting, the physical requirements 
of finite refiectivity and absorption losses by the mirrors. 

Addressing this problem in detail is not only interest- 
ing as a matter of principle, but it is also important for 
a better understanding of the numerous proximity phe- 
nomena arising from quantum fiuctuations of the e.m. 
field described earlier. In many theoretical investigations 
of these phenomena, one deals with conductors that are 
frequently modelled as ideal mirrors. A famous example 
of this is provided by original Casimir's derivation [20| 
of the effect that goes under his name. It is then impor- 
tant to know to what extent conclusions drawn from the 
ideal-metal model can be trusted. Indeed Casimir physics 
offers examples where predictions drawn from the ideal- 
metal model are in contradiction with those derived by 
more realistic modelling of the plates. One such example 
is still much debated as we write, and it is the problem of 
determining the influence of temperature on the magni- 
tude of the Casimir force between two metallic plates in 
vacuum. It turns out the the ideal-metal model predicts 
a thermal force that, for sufficiently large separations be- 
tween the plates, attains a magnitude which is twice the 
one calculated on the basis of realistic dielectric mod- 
els of a conductor, displaying a finite^ though large, dc 
conductivity (for a review of this puzzle, see for example 
Ref. [2]| and References therein). 

In order to shed light on this question, in this pa- 
per we work out a detailed analysis of the full three- 
dimensional commutators for the e.m. field, in the pres- 
ence of dielectric and/or conducting walls with arbitrary 
dispersion and dissipative features. The analysis turns 
out to be considerably more involved than the simple 
one-dimensional model studied in [l3]- Our main result 
is that the canonical commutation relations satisfied by 
free e.m. fields are always valid at all points between 
two macroscopic dielectric or conducting slabs, including 
their surfaces, in full agreement with expectations based 
on the microscopic theory for a system of charged non- 
relativistic particles interacting with the e.m. field. This 
result is consequence of analyticity and fall-off properties 
at large frequencies of reflection coefficients of all real ma- 
terials, as it was correctly conjectured in Ref. [13] ■ We 
also show that such a result is not recovered, however, 
in the case of conductors, if they are modelled as perfect 
mirrors. In this case we find that near the conductors the 
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equal-time commutation relations of the vector potential 
with the electric field have a different form from the free- 
field case. Only at points that are sufficiently far from 
the conductors, the free- fields commutators are recov- 
ered. Our results generalize those obtained by Milonni, 
in the one slab setting, and show that the modified form 
of the commutation relation entailed by perfect-mirror 
b.c. are indeed an artifact of these idealized b.c, not 
shared by real materials. 

The paper is organized as follows. In Section II we 
recall the basic commutation relations satisfied, within 
the microscopic theory, by the e.m. field in vacuum and 
in the presence of charged particles. In Section III we 
briefly review some general results of linear response the- 
ory, as applied to macroscopic quantum electrodynamics, 
and derive formulae for the expectation values of the field 
commutators outside a system of macroscopic bodies, in 
terms of suitable classical Green's functions. In Section 
IV we estimate the Green's functions for a system of one 
or two dielectric and/or conducting slabs in vacuum, and 
in Sec V we use them to calculate the commutation rela- 
tions satisfied by the e.m. field outside the slabs. In Sec 
VI we consider the case of ideal, perfectly reflecting slabs, 
while Section VII contains our conclusions. Finally, three 
Appendices close the paper. 



II. COMMUTATION RELATIONS FOR E.M. 
FIELDS: MICROSCOPIC THEORY 

In this Section we briefly recall well known properties 
of the commutation relations satisfied by e.m. fields, in 
the framework of a microscopic theory of non-relativistic 
matter, where ponderable matter is modelled as a collec- 
tion of non-relativistic charged particles. Here and after- 
wards, we work in Gaussian e.m. units, and we adopt the 
Coulomb gauge. As it is well known the Coulomb gauge 
is very convenient for studying problems where matter is 
non-relativistic, and high-energy processes are neglected, 
for it allows a clear separation of electrostatic and mag- 
netic couplings. In this gauge, quantization is straight- 
forward (see for example the book [23] ) ■ We consider flrst 
the case of free-fields. 



A. Free fields 

In empty space. Maxwell equations imply that the elec- 
tric field is purely transverse 

1 9Aj, 
dt ' 



[A^,{r,t),E^j{v',t)] = -i7Tinc5fJv^r') , (4) 



E 



(1) 



where A^ is the transverse vector potential: 

V • A^ = . (2) 

The fields A_l and Ej^ satisfy the following well known 
equal-time canonical commutation relations: 



where (^^(x) is the transverse delta function [2§|: 



(5) 



(6) 



with k = |k|. 

B. e.m. fields coupled to charged particles 



When charged particles are present, the phase space 
of the total system includes, besides the transverse e.m. 
fields A_L and E_l, the positions x("\ the conjugate mo- 
menta p^"^ and the spins s^") of the particles (labelled 
by the index a). They satisfy the standard (equal-time) 
commutation relations of non-relativistic Quantum Me- 
chanics: 



J") J/3) 



(7) 



(8) 



with all other commutators vanishing. In particular 
x^°'\ pf^ and s|"^ all commute with the transverse e.m. 
fields A_L and Ej^. Finally Aj_ and Ej^ satisfy the 
same equal-time commutation relations holding in empty 
space, Eqs.(l3][5l). 

When charges are present, the electric field E has also 
a longitudinal component Ep: 

E = E||+E_L, (9) 

where E^ is still given by Eq. ([1]), while E|| is equal to 

E„(r,t) = -VC/(r,i), (10) 

where U is the scalar potential: 



U{v,t) 



x(")(t)| ' 



(11) 



with e^"^ the charge of particle a. The scalar potential 
has to be regarded as a function of the particles posi- 
tions, and therefore it is not an independent degree of 
freedom of the system. Since the particle positions com- 
mute among themselves and with the transverse fields, 
Aj^ and Ej^, it follows that 



(3) 



[C/(r,t),C/(r',i)] = 0, 



[U(v,t),A^,{v',t)]=Q, 



(12) 



(13) 



(14) 
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The above Equations imply that the equal-time commu- 
tation relations Eqs.([3][5]) remain valid, irrespective of the 
number and positions of the charged particles, if we re- 
place everywhere the transverse electric field E± j by the 
total electric field Ej : 

[v4^,(r,i),A^,(r',i)] =0, (15) 
[A^,{r,t),E^{r\t)]^-47Tihcd^^{r-r') , (16) 

[E,ir,t),Ej{r',t)]=0. (17) 

The obvious conclusion that can be drawn from these el- 
ementary remarks is that, within the microscopic theory, 
the canonical equal-time commutation relations satisfied 
by the e.m. fields, Eqs. (I3][5]) or, alternatively, Eqs. p^ 
\17\ should be valid always, and therefore they should 
hold, in particular, inside a cavity made of an arbitrary 
material. 



III. COMMUTATION RELATIONS FOR E.M. 
FIELDS: MACROSCOPIC THEORY 

In this Section, we recall a few basic formulas from 
linear-response theory and we discuss the type of probes 
that are needed in order to obtain the commutation rela- 
tions satisfied by the macroscopic e.m. field in the pres- 
ence of dielectrics and conductors. For a review of linear- 
response theory we address the reader to Refs.f?!. 

In linear-response theory, one considers a quantum- 
mechanical system, characterized by a (time- 
independent) Hamiltonian Hq, in a state of thermal 
equilibrium described by the density matrix p 

p ^ e-^" /tr{e-'^") , (18) 

where (3 = l/(fcBT), with ks Boltzmann constant and 
T the temperature. The system is then perturbed by an 
external perturbation of the form: 

Ho.t ^ - fdhY, Qj (r, t) f, (r, t) (19) 

where fj{r,t) are the external classical forces, and 
Qj{r,t) is the dynamical variable of the system conju- 
gate to the force fj{r, t). One may assume, without loss 
of generality, that the equilibrium values of the quanti- 
ties Qj{r,t) all vanish: {Qj{r,t)) = 0. The presence of 
the external forces causes a deviation 6{Qi{r,t)) of the 
expectation values of Qj(r,t) from their equilibrium val- 
ues. If the forces fj{r, t) are sufficiently weak, d{Qi{r, t)) 
can be taken to be linear functionals of the applied forces 
/j(r,t), and one may write: 

6{Q,{r,t)) ^Y. I / dt'4>,,{v,v\t-t')f,{v'X) ■ 

(20) 



The above Equation assumes that the system was in equi- 
librium at t = — oo, and that it reacts to the external 
force in a causal way. The quantities (l)ij{r,r',t — t') 
are called response functions of the system. In princi- 
ple, they can be measured by applying to the system of 
interest suitable external classical probes. 

By a straightforward computation in time-dependent 
perturbation theory one may prove that the response 
functions 4>ij{r,r',t — t') are related to the equihbrium 
(i.e. in the absence of the external forces) expectation 
values of the commutators of the dynamical variables 

Q^(v,t)■. 

0,, {vy,t-t')='-{ [Q, (r, t) , Q, (r', t')] ) e{t -t'), (21) 

where 9{x) is Heaviside step function {6{x) = 1 for x > 
0, 9{x) = for X < 0) and Qi{r,t) is the Heisenberg 
operator: 

Q,(r, t) = e'^^'/'^Q^ir, 0)e-*^°*/^ . (22) 

As it is well known, Eq. (|2ip is the starting point from 
which several general fluctuation-dissipation theorems 
can be derived, that allow to express the (symmetrized) 
correlation functions of the quantities Qi(r, t) in terms of 
the dissipative component of the response functions . 
Since we shall not make use of these theorems in what 
follows, we shall not present them here, and we address 
the interested reader to Refs.fTl for details. 

We wish to exploit Eq. (|2ip to study the commu- 
tation relations satisfied by the macroscopic e.m. field 
at points placed outside a number of dielectric or con- 
ducting bodies. For this purpose, following Agarwal p^ . 
we take the external probes to be a system of classical 
electric and magnetic dipoles, with densities P(r,i) and 
M(r, t) respectively, placed outside the bodies. The ex- 
ternal Hamiltonian iJoxt is then of the form: 

iJext ^-J d3r[p(°'^*)(r,t).E(r,t)+M(^'^*)(r,t).B(r,i)] . 

(23) 

It is convenient for our purposes to have distinct probes 
for the longitudinal and the transverse components of 
the e.m. field. This can be achieved by demanding that 
p(ext) ^^^Y free 

V X P^^'^*) = , . (24) 

If we now express in Eq. (|23p the e.m. field in terms of 
the scalar and vector potentials: 

E = -WU--^ , B = VxA^, (25) 
c ot 

after an integration by parts, and exploiting Eq. (|24[) . 
the external Hamiltonian can be rewritten as: 

i/ext = J d3r[C/(r,i)p(-')(r,i)-iA^(r,t).jf*)(r,t)] , 

(26) 
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(oxt) 



where = -V-P('='^*) andj 

that the current j^'^*'' is transverse: 



Vj 



(oxt) _ 







cVxM(<='^*). Note 



(27) 



We remark once again that the scalar potential U{r,t), 
in the external Hamiltonian Eq. (|26p does not repre- 
sent an independent dynamical variable, and it must be 
regarded as a function of the particle's position, accord- 
ing to Eq.(fTT|). Therefore, in the absence of matter, no 
such term is present in the external Hamiltonian, and the 
scalar potential is zero. 

The response functions are then computed by solving 
the classical macroscopic Maxwell Equations with 
and jIi^'^*'' as external sources: 

V • D = Anp 



(cxt) 



VxE + i^^O 
c at 



V-B = 0, 
V V TT - ^ ^° 4- ,(cxt) 



(28) 



(29) 



(30) 



(31) 



where 



.(cxt) .(cxt) 

Jtot ^ J± 



^p(cxt) 

dt 



The above equations have to be solved subject to the 
usual b.c. of macroscopic electrodynamics, namely (i) 
tangential components of E and H and (ii) normal com- 
ponents of D and B must be continuous across the bodies 
interfaces, which are assumed to have sharp boundaries. 
Coherently with the spirit of a macroscopic approach, the 
dielectrics and the conductors will be described in terms 
of the appropriate electric and magnetic susceptibilities. 
We suppose from now on that the bodies are made of non- 
magnetic (/i — 1)), isotropic and spatially non-dispersive 
materials, characterized by a frequency-dependent elec- 
tric permittivity e(r,a;). We also assume that the bodies 
are homogeneous, in such a way that the permittivity 
e(r, w) is independent of r within the volume occupied 
by each body, with discontinuities occurring only at the 
bodies interfaces. By virtue of homogeneity of the bod- 
ies, and of linearity of the b.c. at the bodies interface, the 
field Equations for the scalar potential U{r, t) are com- 
pletely decoupled from those for the transverse vector 
potential A(r,t). Therefore, we have two independent 
sets of Green's functions: 



U{r,t) 



dt' / d3r'G(r,r',i-t')p^°"*ni-',i') 



Aj.(r,t) - dt' I d3r'G^(r,r',t-t')-jf *^(r',0 , 

(34) 

where Gj^(r, r', t — t') has to be understood as a dyadic 
Green function. 

From the general result of linear-response theory, Eq. 
(I^Tj) . we then obtain the following expressions for the 
two-times expectation values of the commutators of the 
e.m. potentials: 

{[U{r,t),U{r' ,t')]) = thG{r,r' ,t - t') , (35) 
([C/(r,i),Aj.,(r',t')]) =0, (36) 

{[A^r{r, <), A^,{r', t')]) ^ -ih Gx,, (r, r', < - t') , (37) 

where t > t' . For our purposes, it is convenient to split 
the Green's functions, outside the bodies, as sums of an 
empty-space contribution plus a correction arising from 
the material bodies: 

G(r , r' , i - i') = G(°) (v-v',t-t')+ F^"^"*) {v,v',t-t') , 

(38) 

and 



Gi [v,v',t-t')^ G^l^ {r-r',t-t') + Ff {r,r',t-t'). 

(39) 



Here, G*^"^ and Gj^ denote the Green's functions in free 

space, while f and f||"^'**'' describe the effects re- 
sulting from the presence of the bodies. Such a split- 
ting presents the advantage that all singularities are in- 
(32) eluded in the free parts G'"' and G^^\ while the quan- 
tities F^™^'^^ and F^'^*'' are smooth ordinary functions 
of r and r'. The free-field Green's functions have the 
following well-known expressions: 



and 



G 



(0) 



G<-°\r~r') 



d^k 



1 



r — r 



nS{t-t'), 



(40) 



27r2fc 



„ik.(r-r') 



sin[fcc(t-i')] ■ 
(41) 

The factor 5{t — t') in the expression of G'^°-' expresses the 
instantaneous character of the longitudinal electric field 
in the Coulomb gauge. The expressions for the equal- 
time commutators of the e.m. fields are easily derived 
by taking suitable limits of Eqs. (|35ll37p and of their 
time derivatives, for t t'^ . Upon using Eqs. ([55)1 and 
([55)1 , and exploiting the following three relations that are 
obvious consequences of Eq. ([1]): 



1 d 



([C/(r, i), i?^.(r', t')]) = --^([C/(r, t), A^.(r', t')]) , 



([A^,(r,i),i?x,(r',0]) 



1 d 



(42) 



(33) 



^^-([A^.(r,i),A^,(r',0]), 
(43) 
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and 

(44) 

from Eqs. (|35ll37p we obtain: 

{[Uiv,t),U{r',t)])^th A(--')(r,r'), (45) 



([[/(r,i),i?i,(r',t)]) = 0, 



([A^,(r,t),Ax,(r',i)]) = -zhA^'^iry), (48) 



([A^,(r,i),i?j.,(r',i)]) 



and 



i h 



r2 ^ij 



where we defined 



yl("^^')(r,r') = lim ^('"''^Vr, r', t - t') ■ 



(r,r')^^Um f^:;-)(r,r',t-t') 



i3i'^f)(r,r')^^Um ^^(r,r',t-0, 



and 



d^^'\r,r')= lim 



g2^(mat) 



t^t'+ dtdt 



(46) 



(47) 



(49) 



C^fM^ (50) 



(51) 



(52) 



(53) 



(54) 



By comparing Eqs. (|45ll50p with Eqs.(1311Sl) and Eqs. (fT^ 
[T4)) . we see that outside the bodies the free- field canonical 
commutation relations are recovered provided that the 
quantities A^'^^'\ A^J'^f, B^^f and are zero. We 

shall prove below that this indeed the result of 

analyticity and fall off properties at large frequencies of 
the reflection coefficients of all real materials. 

Before we turn to detailed computations, we present 
below the field Equations satisfied by G and G_l. They 
are conveniently expressed in terms of the (one-sided) 
Fourier transforms of the Green's functions, defined as: 



G(r,r',w)= / diG'(r,r',i)e'"'* , (55) 
Jo 



G_L(r,r',tj)= / dtGi(r,r',i)e*'^* . (56) 
Jo 

From Maxwell Equations we then obtain: 

V- [e(r,w) VG] -47r5(r-r') , (57) 



(A-f e(r,c^)coVc')Gi(r,r',c^) = -47r <5i(r - r') , (58) 

where 5±{v—v') is the transverse delta-function dyad, Eq. 
([S]). These Equations must be solved with the appropri- 
ate b.c. at the bodies interfaces, and must be subject to 
the conditions required for a retarded Green's function 
[2^. For later use, it is useful to recall the main prop- 
erties enjoyed by the Green's functions [l^]. First of all, 
they satisfy the following reciprocity relations: 



and 



G(r,r',c.) = G(r',r,c.), 



Gii-,-(r, r',w) = Gj_ji(r',r,Lj) , 



(59) 



(60) 



that are a consequence of microscopic reversibility. The 
next set of properties express reality features of the 
Green' functions, and are a direct consequence of real- 
ity of the external sources 



and 



G*(r,r',c.) = G(r,r',-c.) 



G'Lj(r,r',t^) = G±y(r,r',-t. 



(61) 



(62) 



The next set of properties is a consequence of the fact 
that the permittivity e{u)) of any causal medium is an an- 
alytic function of the frequency w in the upper complex 
half- plane C"*" [2^ (see also Appendix B)). This implies 
that the Green's functions G(r,r',Ci;) and Gx(r, r',w) 
are also analytic in C^, as it must be case for a re- 
tarded response function. In C"*" they satisfy the con- 
ditions G*(r, r',w) = G(r,r',-w*) and G*^^j{v,Y' ,w) = 

Gx,;j(r, r', — w*), that generalize the reality conditions 
Eq. (pT|) and Eq. respectively. These more gen- 

eral properties imply that the Green's functions are real 
along the imaginary frequency axis: 



and 



G(r,r',zO=G*(r,r',zC), 



Gi.,(r,r',zO=Gl„-(r,r',iO 



(63) 



(64) 



It is finally useful to write down the inversion formulas 
expressing the Green's functions, in the time domain, in 
terms of their Fourier transforms. They are 



G(r,r',i-i') = / dwG{v,v',w)e 



-iw(t — t') 



(65) 
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G±{r,r',t-t') = j rfwGi(r,r','u;)e-™(*-*') 



(66) 



where T is any contour in C+ that can be obtained by 
smoothly deforming the real frequency axis, keeping fixed 
the end-points at infinity. Analyticity of the Green's 
functions in ensures that the integrals on the r.h.s. 
are independent of the chosen contour F. 

In the next two sections we shall compute the Green 
functions at points outside a single dielectric slab, and 
between two plane parallel slabs. 



IV. GREEN'S FUNCTIONS OUTSIDE 
DIELECTRICS AND CONDUCTORS 

In this Section we evaluate the e.m. Green's functions 
outside dielectric and/or conducting slabs. 

In the next two subsections we shall separately consider 
the cases of one slab in vacuum, and two plane-parallel 
slabs separated by an empty gap. We choose our carte- 
sian coordinate system such that the z-axis is perpendic- 
ular to the slabs. Translational invariance of the system 
in the [x, y) plane implies that the quantities F^™*^*) and 
-p(mat) functions only of z, z' and (r^ — r'^), where 
we denote by xj_ the projection of the vector x onto the 
(x, y) plane. The computation is facilitated if we express 

G'"' and G*^'' in a form that is adapted to the symme- 
tries of our problem. Consider first the free scalar Green's 
function G^") : 



G(°)(r-r') = 



1 



(67) 



We note that G'-^-' is independent of the complex fre- 
quency w, as it must be because of the instantaneous 
character of the longitudinal electric field in the Coulomb 
gauge. For our purposes, the convenient form of is 
the following well-known Weyl representation: 



(5(0) ^ 



27rfci 



2k_L ■ (r_L — rj_ ) — \ ^^^' \ 



(68) 



that can be easily obtained by integrating over ky, the 
standard plane- wave decomposition of G'-^-' . The above 
expression for can also be written as 



G(o) 



(69) 



2ti k± 

where we define k*-^' = ±ikj_ z, and the upper (lower) 
sign is for z > z' (z < z'). Consider now the familiar 

representation of G*^"*: 
d^k 1 



27r2 fc2 _ ^2 I 



'ik-(r— r ) 



(70) 



where fco = w/c. In Appendix A, we show that G^^ can 
be decomposed as the sum of two dyads: 



Gf = U(") + V(") 



(71) 



Here, U^^^ denotes the tensor of components 




where kz = ^JkQ — k'^ (the square root is defined such 
that lm{kz) > 0), = z x k^, k^*' = kj^ ± k^z and 
= kj_z ^ kzii±. As to Vi\^\ it can be written as 



(0) 



1 92^(0) 



where ^'^"^ is the function 



(0) 



2tt k. 



^ g4k(±'-(r-r') 



(73) 



(74) 



In both Eqs. ([7^ and ([71)) the upper (lower) sign is 
for z > z' (z < z'). It is useful to provide a simple 
intuitive interpretation for the above Green's functions 
that will be useful later when we consider the influence 
of a material slab. Consider first the expression for G'^"-' 
given in Eq. (j69p : we can interpret is as consisting of 
a superposition of instantaneous scalar waves originat- 
ing from point r', that propagate to the right (left) with 
wave- vector k*-^) (k^"^). Consider now our expression 

for Gf \ Eq.dni). Its first contribution U^"), Eq. ([7^. 
can be physically interpreted as a superposition of e.m. 
waves with TE and TM polarization corresponding, re- 
spectively, to the first and second term between the round 
brackets in Eq. (j72p . These waves originate from point 
r' and propagate to the right (left) with wave-vector k(+^ 
(k^^)). We note that for feg > k±_ these modes represent 
propagating waves, while for fco < k± they are evanescent 
waves, that decay exponentially as we move away from 
z'. The second contribution to G^-* , V^'^^ , can instead be 
interpreted as representing scalar waves that propagate 
instantaneously from point r' in the right (left) direction, 
with wave- vector k(+^ (k*^^)). 

We are now ready to compute j^(bodios) -p(bodics) 
We consider first the one-slab case. 



A. The case of one slab 

In this Section we compute the Green's functions out- 
side a single dielectric or conducting slab, occupying the 
half-space 2 < 0. Following the remarks of the previ- 
ous Section, outside the slab and on its surface, i.e. for 
z,z' >0 we define: 

G("^")(r,r',u.) = G'^°\r-r')+ F'^'"^^^\r^-r'^.z.z',w) , 

(75) 

and 

G^""")(r,r',u;) = G^°\v-v' ,w)+¥;^^''\r^-r'^, z, z' ,w) 



(76) 
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Fixine; once and for all z' > 0, we make for F^™'^ the 
following ansatz: 



e[w)k^ - q 
£{w)k^ + q 



(85) 



F 



(wall) _ 



27rfc 



r(w) e-- ir±-rl)-kAz+z') . > . 



(77) 

For z < 0, the complete Green's function is taken to be 
of the form: 



(5 (wall) 



271 kA 



t{w) e^k_L-(ri-r^)-fci(2 ^ z<0. 

(78) 

Both ansatz ensure appropriate fall off for |z| ^ cxd. It 
is easy to verify that the above ansatz satisfy the b.c. at 
z = 0, provided that we take: 



r{w) 



e{w) - 1 
e{w) + 1 



and 



t{w) = 1 — r{w) 



(79) 



(80) 



The chosen forms of F^"^"), for z > 0, and G^™'^), for 
z < 0, have a simple physical interpretation, that will 
be be useful later when we shall consider the more elabo- 
rate case of two slabs. In empty space, the source p(r', it;) 
generates "instantaneous" scalar waves of (complex) fre- 
quency w originating at r' and propagating in the right 
direction (i.e. towards larger z) with (complex) wave- 
vector k'+'', and in the left direction with wave- vector 
k*^"^. When a wall is present, the left-moving waves hit 
the wall and then we have a reflected wave with ampli- 
tude f{w), and a transmitted wave of amplitude t(w). 

We can now evaluate F^'^'^-'. In a way analogous to 
Eq. ([7T|) . we decompose it as: 



-p,(wall) _ •jj(wall) Y-(wall) 



(81) 



Inside the slab, for the full Green's function we set in- 
stead: 



j=;(wall) ^ ^(in) ^ ^(i„) ^ 



(82) 



Linearity of the boundary-value problem permits to de- 
termine separately U(™'^) and V^™''^ The physical pic- 
ture of U^°) as a superposition of TE and TM waves, 
suggests at once the following ansatz for U^^^''); 



(wall) 



2 



— (e±ie±j r^*^ {w, k±) 



where q = e{w)kQ — k^. A somewhat lengthy solution 
of the boundary- value problem indeed confirms the above 
intuitive form of J/^'^''''^. Consider now V^*^"). Eq. ^ 
suggests that we set: 



(wall) 



1 g2^j;(wall) 



while inside the slab (i.e. for z < 0) we set: 



(in) _ 



1 52^J,(in) 

k^ dxidx'j 



(86) 



(87) 



It can be seen that the appropriate boundary dielectric 
conditions at z = are satisfied, provided that the func- 
tions ^(^^^^'> and fulfill there the following b.c: 



\z=0 



^(wall)^ 



2=0 



e(w)*('")'| 



(^(0)' + |r(wall)'^ 



(89) 



where a prime denotes a derivative with respect to z. 
One then finds: 



(^(wall) ^ _ 



27r kj_ 



-{w) e 



ik'+'-r-ik'-'-r' 



(90) 



where r(a;) is the reflection coefficient in Eq. ([79]) . We 
note that the expression of ^pt™'') coincides with that of 
^(waii)^ We remark that F^*^") and f|l^''"' are analytic 
functions of the frequency w in the upper complex plane 
C~^, as a result of analyticity in C+ of the reflection coeffi- 
cients fiw), {w) (see Appendix B). Moreover, we note 
that has no singularities along the real-frequency 

axis, as it can be easily checked from Eq. ((77|) . if one 
considers that the refiection coefficient f{w) is finite in 

(see Appendix B). As to F^'^'^-', it only has an in- 
tegrable singularity at kz =0. The presence of singular 
factors proportional to k^^ in the expressions of U^J'^^^'^ 



and v}""^^^^ (see Eq. 



and Eq. 



does not cause 



any further singularities at w = 0, for it can be verified 
that these singular terms cancel each other upon taking 



the sum of t/^,™") andyA"'''^ as we now show. Indeed, 
upon collecting in Eq. ([53)1 and Eq. ([55)1 the terms that 
are singular at w = 0, we obtain: 



^(wall) 



ko J 

where r*^") (w, kj_) and r'-^^(w,kj_) are the familiar Fres- 
nel reflections coefficients for TE and TM waves, respec- 
tively: 



r^'\w,k±) 



kz -q 
kz + q ' 



(84) 



lim 



(r,r')) 



/ 27rfc I 



,,2 t '1 



(91) 



where we made use of the following relations 

^(±) =^ik(±)+0(w2) , (92) 
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k, = tk^ + 0{w^), (94) 

to substitute everywhere S.'^^K k'^'^' and by ^zk'^^ 
ic^^^ and k± respectively. Now in Appendix B it is shown 
that both, for dielectrics and conductors, the difference 
j-(p) — f approaches zero as w^: 

r(p){w) -f{w) ^0(w^) . (95) 

Therefore, the ratio (r^^*' — f)/w^ is finite as ui tends to 
zero, showing that -Pjjj'"'' is regular at w = 0. 

B. The case of two plane-parallel slabs 

In this Section we calculate the Green's functions for 
the case of a cavity constituted by two non-magnetic ho- 
mogeneous, isotropic and spatially non-dispersive plane- 
parallel slabs separated by vacuum. We assume that the 
slabs can be characterized by the respective electric per- 
mittivities, ei{w) and 62 (w). We choose our cartesian 
coordinate system in such a way that slab one occupies 
the region — 00 < z < 0, while slab two occupies the re- 
gion d < z < 00, d being the separation between the two 
slabs. The formulae derived in the preceding Section, for 
the one slab case, can be easily generalized to the two 
slabs setting, on the basis of the intuitive physical pic- 
ture of the free Green's functions as consisting of left and 
right moving waves originating from r' . 

Let us consider first the scalar Green's function G. 
Analogously to what we did in the previous Section, in- 



r7(cav) _ [ c^^ki. 



1 



where r^^"'', a = s,p are the Fresnel reflection coef- 
fcients of slab i for polarization a, and Aa = 1 — 
r^"V^"^ exp(2ifc^d). For V^'^^^) we obtain: 



side the cavity (i.e. for < z, z' < d) we set: 

&'"'^\t,t' ,w) ^&^\t-v')+ F'^''''^\r^-v'^,z,z',w) . 

(96) 

The expression that one finds for F'^'^'^^'^ is analogous to 
^^(waii)^ but of course one must take account now of the 
possibility of multiple reflections off the two slabs. This 
is easily done, by inserting for each reflection by slab i 
the appropriate reflection coefficients fi{w), that has an 
expression analogous to Eq. (|79p (with ei{w) in the place 
of £{w)). Moreover, a factor e^'^^^'^ must be included for 
each round-way trip from one slab to the other and back. 
One obtains: 




(97) 

where A — 1 — ri(w) r2{w) exp(— 2fc_L d). 
For the transverse Green's function, we set: 



G^'"")(r,r',w;) = Gf {r-r' .w)+f'^l'-^\v ^-v'^, z, z' ,w) , 

(98) 

with 

where the symbols have the obvious meaning, analo- 
gously to previous Section. The same arguments that 
led us to write Eq. ([57]) now give: 



(100) 

I 

where 



1 

~As 



''l ik'+'-r-ik'^'-r' , ^2 ik<"' ■r-ik<+' ■r'+2i/c^(i 

A A 



— -l \ (^^(+)^(+)e»k<+'-(r-r') _^^(-)^(-)g»k<-).(r-r')'j , A+) A-) ^1 ' ^»k<+) ■r-»k'-)r' 



(-)c(+) ^2" „ik<-)-r-zk'+>r'+2tfc^rf 



(101) ^ J ^^^-^ [^1 - 1^ |^gik<+) (r-r') ^ g»k<-).(r-r')^ 
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ik<+'-r-ik<"'- 



+ r2 e 



ik<-'T-ik<+'T'-2fe I d 



(102) 

Again we find, as in one slab case, that the expression 
of \l/(cav) coincides with that of F^'^'^^\ The same consid- 
erations used in the one-slab case can be now repeated 
for and F^'*^-' to show that both quantities are 

analytic in C+, and have a finite limit for vanishing w. 



V. COMMUTATION RELATIONS FOR THE EM 
FIELDS INSIDE A CAVITY 

In this Section we compute the quantities ^f™*^'), 
^(inat) ^ ^(™^*) a,nd C'f^j^^ for the two slab setting consid- 
ered in the previous Section. The corresponding quanti- 



ties shall be denoted by A^^j\ B^l'^^^ and C 



±ij 



±ij 



i(cav) 



respectively. We shall see that they all vanish, as a con- 
sequence of the analyticity and fall-off properties at large 
frequencies of the reflection coefficients of all real mate- 
rials. As seen in Sec. Ill, vanishing of these quantities 
entails that the e.m. field satisfies free-field commutation 
relations in the empty region between the slabs. 

Consider first the quantity A^'^'^^\r,r'). From its defi- 
nition Eq. ([5T|) it follows that can be expressed in 
terms of as: 



'(r,r')= lim / ^F^'^^^) (r, r', lu) e"™ ^ , (103) 

T^0+ Jp ZTT 



where is given in Eq. ([57)1 . In Appendix B it 

is shown that the reflection coefficient f(w) of any real 
material vanishes like w~'^ for large values of \w\ and 
this implies, as can be seen by inspection of Eq. (|97p . 
that F^^''^) approaches zero like w^"^. Therefore is 
absolutely integrable, and then in Eq. (|103p we can take 
the r-limit inside the integral. After we do it we obtain: 



Aic^^) (^r,r') = f — #('=^^)(r,r',«;). 
Jr 27r 



(104) 



The w fall-off rate of F'^'^'^'^^ at infinity now permits to 
close the integration contour in Eq. (|104[) in the upper 
complex w-planc C+, and then analyticity of F^'^^^) in 
implies at once that the integral is zero. Therefore 
we conclude 



yl('=^^)(r,r') =0 



(105) 



We turn now to the quantity (r, r'). In view of 



its definition Eq. ([5^ we have 

(r, r') = ^Inn ^Fj^^ (r, r', w) e"™ ^ , (106) 

and then to prove that it vanishes, we need consider the 
fall-off properties of F^fi^^- According to Eq. ([99]), it is 

the sum of two terms: fI'T' = U^^''"'^ + v'^T^ . As to 



^iTj''' we see, by inspection of Eqs. ((TOT|) and (fT02)) . that 
the fall-off rate of f{w) implies that 1^^'^^^'' falls-off like 

w^*. Consider now u'^l^j'^ ■ We note first that, because of 
the kz factors in the exponentials, all terms in the r.h.s. of 
Eq. ()100|) decay exponentially fast as w goes to infinity 
in along any direction not parallel to the real axis. 
Along the real axis, since Fresnel refiection coefficients of 
all real materials decay like (see Appendix B), u'f^^j'^ 
decays at least as fast as (in fact a more careful 

analysis carried out in Appendix C shows that the rate 
of decay is actually w"**). Therefore, F^^^^ decays in all 
directions in at least like and then, by following 
exactly the same reasoning used in the case of A^'^^'''^ (r, r') 
we can prove that 



A^>{ry)=0 



(107) 



We remark that the above equation holds also when ei- 
ther r or r' or both belong to the slabs surfaces. Consider 
now the quantity B'f^'^J^ (r, r'). Recalling its definition Eq. 
(|53| . we have 



B 



(cav) I 



->0+ 



dw 



(r, r') = i lim / ^ w F\^^'"\r , r' , w) e 



27r 



±ij 



(108) 

Thanks to the w""^ fall-off rate of F^f^^\'w), the extra 
power of w does not spoil convergence of the w-integral on 
the r.h.s. of Eq. (jlOSp . and therefore the same arguments 
used to prove that A^^'^J^ is zero can be used to obtain 



j-){ca,v) , /\ 



0. 



(109) 



Finally, we consider the quantity c'f^j\r,r'). For this 
we have 



C 



(cav) / 
±ij 



(r,r')= lim / ^ F|<:f (r, r', ^ e 



+0+ 



27r 



(110) 

Proving that C^^^J^ vanishes requires much more labor, 
because of the two extra powers of w in the integrand 
on the r.h.s. of Eq. (|110p . We relegate the proof in 

Appendix C, where we show that the decay rate of Fj^^^' 
is actually which is sufficiently fast to imply: 







(111) 



.(cav) n(cav) 



Having proved that the quantities A'^'^ 

and C^f!^^'^ vanish, we then reach the important con- 
clusion that in the empty space between two dielectric 
and/or conducting slabs the e.m. fields satisfy free-field 
equal-time commutation relations, Eq. (|12j|17p . This re- 
sult is consistent with what was expected on the basis 
of the microscopic theory, for a system of non-relativistic 
charged particles interacting with the e.m. field, as we 
have seen in Sec. II. We remark that no singularities are 
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encountered as r and r' approach the slabs surfaces, and 
therefore the canonical form of the free-space commuta- 
tors also holds on the surfaces of the slabs. It is important 
to realize that these results are intimately tied to analyt- 
icity and fall-off properties of the reflection coefficients of 
real materials. 



VI. COMMUTATION RELATIONS OUTSIDE 
IDEAL CONDUCTORS 

In this Section we investigate the commutation rela- 
tions satisfied by the e.m. fields outside ideal conductors. 
Ideal conductors are characterized by the fact that they 
have constant reflection coefficients. Indeed, by taking 
the limit e ^ oo in Eqs. §^ and (gS]), we find 

that for an ideal conductor f and r^*'' are one, and 
is minus one at all frequencies. Obviously, constant re- 
flection coefficients are analytic in C~^, and therefore the 
main difference between ideal conductors and real ones is 
that the reflection coefficients of the former do not van- 
ish in the limit of large frequencies. We shall see below 
that this feature entails that the e.m. field outside the 
conductors, and on their surfaces, fail to satisfy free-field 
canonical equal-time commutation relations. 

In order to determine the commutation relations sat- 
isfied by the e.m. field we consider again Eqs. (|45ll50p 
that remain valid also for ideal conductors. All that we 
have to do then is to evaluate the quantities on the r.h.s. 
of these equations, using the values of the reflection coef- 
ficients pertaining to ideal conductors. We consider first 
the simpler case of a single conducting slab. 

We start by evaluating the quantity jT'CidwaU)^ where 
the superscript (id wall) stands for a slab made of an ideal 
metal. From Eq. ^ we note that for f{w) = 1, p^''^'"''^^^ 
becomes independent of the frequency and, upon taking 
the inverse time-Fourier transform, one easily finds that 
^(idwaii) jg proportional to S{t - t'). Then F^''^^'^") is 
zero for all t > t' and therefore from Eq. ((35)) we have 

([[/(r, t), U{r', t')]) = (ideal conductors). (112) 

Upon taking account also of Eq. ([M]) we see that outside 
an ideal conductor, all two-times commutators involving 
the scalar potential U have vanishing expectation values, 
and this implies 

U{r, t) =0 (ideal conductors). 

Therefore, outside an ideal conducting slab the longitu- 
dinal electric field is zero. We evaluate now the quantity 
^I'ij^^''^- Upon using the identity 



A+)A-) 



-XkS^kSjk + " J , (113) 

Kg 



where Ai = A2 = —A3 = 1, one finds that, for r = r^'^^ 



1 and r^-^^ = - 

^(id wall) 



F 



1, Eqs.dSI]), ([831) and (|90D lead to: 



27r 



Y I ^k^^k^Jk p 



"0 



^ ^ i _^-k^{z+z') 



fr2 



(114) 

By a similar computation as the one described in Ap- 
pendix A, it is possible to verify that the r.h.s. of the 
above equation can also be written in the following form: 



^(id wall) 



1 



X e 



i\k.±-(r±-r'^)+k3(z+z')] 



(115) 



where k^'') 



fcaz. From this expression we see 



that in the case of an ideal wall p^^^^"^^^ decays for large 
frequencies only like w~'^ , and not like as we found 
in the case of a slab made of a real material. This fall- 
off rate is sufficient to prove, by the same steps used 
in the previous Section, that Eq. (I107|) remains valid. 
Therefore we find that also in the case of an ideal slab 
the equal-time commutators for the vector potential have 
the canonical form Eq. ^ at all points outside the slab, 
including its surface. The w"^ fall-off rate is not sufficient 
however to ensure validity of Eq. (|109p . and we show now 
that for an ideal conductor Eq. (jl09p indeed fails to be 
true. To see this we take the inverse Fourier transform of 
Eq. (|115p , as defined in Eq. . The frequency integral, 
for t > t' can be easily evaluated by closing the contour F 
in the lower complex plane (which is possible now because 
the r.h.s. of Eq. (|115p is analytic also there), and by 
noticing that the integrand has poles only at feg = ±A:. 
We get: 



F 



{id wall) 



27r2fc 



^k^ikSjk 



fc2 



X g^[k_L■(r_L-rl)+fc3(^+^')] sm[ck{t - t')] . (116) 

Then, from Eq. (|49)) we obtain: 

{[A^^ir,t),E^,{r',t)]) 

= -Amhc[S^^ir-r')- <5|f """^ (rj.-rl, z+z')] , (117) 



where we defined 



^kSikSjk 
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J[k±-{r±-r'^) + k3{z+z')] 



(118) 



We note that Sl"^ is a smooth function for z + z' > 
approaching zero for large z and z', but it is singular 
when both z and z' belong to the slab surface (i.e. for 
z = z' = Q). In particular, for i = j = 1, Eq. ()117p gives 



xe 



ik3{2d+z-z') 



= -47ri?ic[(5i^i(r ~r') - (5i^i(ri - r^, z + z')] , (119) 

in agreement with the finding of Ref. [15[. By using 
Eq. (|44|) . and Eq. (|116p . it is easy to verify that the 
canonical commutation relations for the components of 
the transverse electric field, Eq. remain valid. 

We turn now to the more elaborate case of two plane- 
parallel ideal slabs. We shall be brief here, the analysis 
being similar to the one slab case. First we note that, 
similarly to F^'*^™''), also the quantity f'(idcav) becomes 
independent of the frequency when perfectly reflecting 
slabs are considered, as it is easily seen from Eq. ([97|) . 
Therefore, j^C'dcav) jg proportional to 6{t — t'), and again 
we conclude the the scalar potential can be taken to be 
zero outside the slabs. We consider now the transverse 
Green's function. A somewhat lengthy, but straightfor- 
ward computation analogous to the one done for the one- 
slab case, gives the following expression for the quantity 



F 



(id cav) 



p(id cav) _ _ 



d^k 1 

27r2 k yl(id) 



ky>k 



fc2 



^ik3{2d—z+z') 



\ ^ik3{2d+z-z') 



^kSikSjk 



fc2 



■3 1 ^ik3{2d-z-z') 



X gi^^i'^^-'^l) sm[ck{t-t')] 



(120) 



where ^^"^^ = 1 — exp(2i/c3(i). By using this Equation, 
and recalling Eqs. (|48|) and ([50|) . we easily see that the 
equal-time commutators for the vector potential on one 
hand and for the transverse electric field on the other, 
both vanish inside the cavity and on the slabs surfaces, 
in agreement with the free-field case, Eq. ^ and Eq. 
([5]). On the other hand, from Eq. we get: 



-47rznc[4(r-r')- ^^^f ^^^^ (r^ - r^, z, z')] , (121) 



where 



s-(id . 



d^k 1 
(27r)3 yl(id) 



kikj 



+ >^kSik5jk 



k^'^H- 
2^1 22 

fc2 



^lk3{2d-Z-z') 



(122) 

We note that the first and the second terms between 
the square brackets on the r.h.s. of Eq. p22p represent 
smooth functions of z and z' at all points between the 
slabs, including their surfaces, while the third and fourth 
terms are singular, respectively, on the surface of slab one 
(i.e. for z = z' = 0) and slab two (i.e. for z = z' ^ d). 
Moreover, we observe that in the limit of large separa- 
tions d, and for fixed z and z', the phase factors involving 
d in the first, second and fourth terms between the square 
brackets on the r.h.s. of Eq. (|122p . oscillate more and 
more rapidly, and so suppress the corresponding terms. 



In this limit S, 



(id cav) 



tends to S, 



(id wall) 



and then Eq. ([^1]) 



reproduces Eq. (|117p . 

From the above analysis, we see that while all other 
commutators have the free-field form, the presence of the 
extra terms S^"^ ""^^^"^ and 5^"^ '^'^^^ on the r.h.s. of Eqs. 
(|117p and (|12ip respectively, implies that perfect-mirror 
b.c. lead to equal-time commutation relations for the 
vector potential and the electric field of a different form 
from the free- field ones, Eq. ([4]). In the one-slab case, 
Eqs. prz)) and pT5)) show that the free-field form of the 

commutators is recovered only when the quantity S^^^'^ '^''^^ 
can be neglected, and this occurs at points z and z' that 
are far from the slab. In the cavity setting, Eqs. (|12ip 
and (|122p show that free-field commutation relations are 
recovered only provided that the quantities 2c? -I- z — z', 
2d — z + z', z + z' and 2d — z — z' are simultaneously 
large. This is only possible for large cavities, and for 
points z and z' far form both conductors. This is in 
contrast with what was found in the previous Section, 
where we proved that in the case of real materials free- 
fields equal-time canonical commutation relations retain 
their validity everywhere between the slabs, including on 
their surfaces. 



VII. CONCLUDING REMARKS 

In this paper we have determined the commutation 
relations satisfied by the quantized e.m. field outside 
one or two plane-parallel dielectric and/or conducting 
slabs in vacuum, assuming that the slabs are made of 
isotropic and homogeneous, spatially non-dispersive ma- 
terials, with arbitrary frequency-dependent dispersion 
and absorption. Using a general form of macroscopic 
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quantum electrodynamics, we have found that at all 
points between the slabs, including on their surfaces, the 
e.m. field satisfies canonical commutation relations of the 
same form as in empty space, in full agreement with the 
microscopic theory. This result is a general consequence 
of analyticity and fall-off properties at large frequencies 
satisfied by the reflection coefficients of all real materials. 



We have also shown that free-field equal-time commu- 
tation relations do not obtain outside one or two conduct- 
ing slabs, if the latter are modelled as perfect mirrors, 
because of extra terms that appear in the commutator 
of the vector potential with the electric field. Free-field 
commutators are only recovered at points that are suffi- 
ciently far from the mirror. In the one-slab setting, our 
findings coincide with those obtained by Milonni [15| in 
his investigation on the Casimir effect. Since no such de- 
viation from free-field commutators is found in the case 
of real materials, we draw the conclusion that the mod- 
ified form of the field commutation relations implied by 
perfect-mirror b.c. is an artifact of this idealized model, 
as it was conjectured to be the case in Ref.^ff] on the ba- 
sis of a simplified one-dimensional form of QED. Even if 
the commutator of the vector potential with the electric 
field, being a gauge dependent quantity, is not a phys- 
ically observable quantity, failure of perfect mirror b.c. 
to reproduce the correct free-field form of the equal-time 
commutation relations near the surfaces of the conduc- 
tors, indicates that a certain amount of caution should be 
used when these idealized b.c. are used in investigations 
of proximity phenomena originating from the quantized 
e.m. field in the presence of conductors. 



Before closing, we would like to comment on possi- 
ble generalization of the results derived in this paper to 
other materials, including magnetic materials, and non 
isotropic or spatially dispersive media. Consideration 
of isotropic magnetic materials offers no difficulties, be- 
cause it just requires substituting in our formulae the 
well-known expression for Fresnel reflection coefficients, 
for a medium with magnetic permeability fi. On the 
other hand, it is known today that the general formulae 
Eq. (|35ll37p expressing the expectation values of the field 
commutators in terms of their classical Green's functions, 
are valid for arbitrary media [2?! . |28 l| , and therefore one 
can use them also in the case of anisotropic and/or spa- 
tially dispersive media provided only that one is able to 
determine the reflection coefficients for a slab made of 
these materials. Since reflection coefficients of all media 
are analytic functions of the complex frequency in the up- 
per complex plane, and fall off to zero at large frequencies 
[201, it is therefore expected that (free-space) canonical 
commutation relations remain valid also for these more 
general materials. 



APPENDIX A: WEYL REPRESENTATION OF 
THE TRANSVERSE GREEN'S FUNCTION IN 
EMPTY SPACE 

As it is well known, the time Fourier-transform of the 
Green's function for the transverse e.m. field in free space 
is given by the formula: 

^^'J " J ^2 _ ^2 [^^^J fc2 j e . (AiJ 

An expression analogous to Weyl's representation of the 
scalar Green's function, Eq. (|68p . can be obtained by 
performing the integral over in Ea. (|Al[) . The integral 
can be done easily by suitably closing the contour 
of integration in the complex k^ plane (for z > z' one 
closes the contour in upper half-plane, for z < z' in the 
lower plane). The integral then receives contributions 
from poles in the integrand of Eq. ljAip . which arise from 
two sources. The first one is the factor involving the 
inverse of fc^ — which gives rise to poles at k^ = ikz, 
where k^ — y^k^ — k'j_ (the square root is defined such 
that Im(fcz) > 0). Importantly, the second source of 
poles is the gauge fixing term inside the round brackets, 
proportional to the inverse of fc^, which has poles at ^3 = 
zLik±. We correspondingly split G^"* as the sum of two 
terms: 

= U(") + V(o) , (A2) 

where U*^°^ accounts for the former set of poles, and V*^°) 
for the latter. We find: 

Iji") -± f ^ (s ^'^^^f e^k<±).(r-r') 

~2j kz ki j"" 

(A3) 

where k^^' = ± fc^z, and 

where k^^'^^ = k^ ± ikj^z, and in both Eqs. (|A3p and 
(IA4p the upper (lower) sign is for z > z' {z < z'). It 
is now convenient to further transform the expression of 
^ij ^ by considering the following decomposition of the 
identity 6ij: 




where = z x k^, and ^ = fc^z =p kz]i.±. Upon re- 
placing Sij, in Eq. (|A3p . by the r.h.s. of Eq. (|A5p . we 
obtain: 

~2nJ kz [^^^^^^+ kl ) 

(A6) 
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As for Fj^-"' , we note that it represents a pure scalar con- 
tribution, for it can be written as: 



~ (0) ^ 1 



kl dx.dx'^ 



where 



^ = 



27rfc 



^_ gik(±'-(r-r') 



(A7) 



(A8) 



Upon combining Eq. (jA6P and Eq. (|A7p . we obtain the 
foUowing final expression for G 



A±)A±)\ 

'^i '^i I „ik(±'-(r-r') 



1 d^^ 

+ To 



/cq dxidx'j 



(A9) 



It is useful to consider also the behavior of the reflec- 
tion coefficient in the limit of zero frequency. In the case 
of a dielectric, e{w) approaches a positive constant eg > 1 
at zero frequency, and therefore for the reflection coeffi- 
cients f(0), M{0) and A'^O) we find: 



f(0) 



= r(P)(0) = 



£0-1 
eo + 1 



< 1 , (dielectrics) (B3) 



r('')(0) = (dielectrics). 



(B4) 



Consider now conductors. At sufficiently low frequency, 
the permittivity of a conductor is of the form 

e{w) — A-Ki — , (conductors) (B5) 
w 

where ctq is ohmic conductivity. Then for the reflection 
coefficients at zero frequency we find: 

f(0) = r(P)(0) = 1 , (conductors) (B6) 



APPENDIX B: PROPERTIES OF THE 
REFLECTION COEFFICIENTS 

In this Appendix, we briefly review some general prop- 
erties of the reflection coefficients of dielectrics and con- 
ductors, that are important for the present paper. They 
are a direct consequence of the general properties of the 
electric permittivity e(w). As it is well known ^2^, the 
permittivity of a causal medium is an analytic func- 
tion of the complex frequency w in the upper complex 
plane . Moreover, its imaginary part e" is never 
zero in C^, except along the positive imaginary axis 
(w = ,^ > 0), where e is positive and monotonically 
decreasing. For large (complex) frequencies w, the elec- 
tric permittivities of all materials approach one, and have 
an asymptotic expansion of the form [25l |: 



r('''(0) — (conductors). 



(B7) 



It is also useful to estimate the behavior of the difference 
between Ap\w) and r(w)), as w approaches zero. In the 
case of dielectrics, one finds: 

while, in the case of conductors, we find: 
iw^ / 3 ttctq 



+ 0(w*) . (B9) 



We see that in both cases, the difference r^^^ — r ap- 
proaches zero as w'^. 



e{w) 



1 



A B 



(Bl) 



where A and B are real positive constants characteristic 
of the material. 

As a consequence of the above analyticity properties of 
£(?«), the reflection coefficients f, r'^' and given by 
Eas. (|79| . ((85| and ((84)) respectively, are analytic func- 
tions in C"'", and they are real and positive along the 
positive imaginary axis j26| . The asymptotic behav- 
ior of e{w) implies that for large frequencies f(w) and 
r^"-*(w, kj_), a = s,p have asymptotic expansions of the 
form: 

CD / , C(") D(") 

f{w) ~ , /"^(w,/cj_) ~ — + , 

(B2) 

where C, D, C^"' and L>(") are real numbers character- 
istic of the material. We remark that C*^"^ and are 
independent of k±. 



APPENDIX C: PROOF THAT C^^^j'^ = 



In this Appendix we prove that the quantity C^^^J^ 
deflned in Eq. piOp is zero. 

Upon recalling Eq. (|98p. we first split C^f^^^ as 

cl7)(r,r') = A,(r,r') + i?.,(r,r'), (CI) 



where 



and 



g2jj(cav) 



^2y-(cav) 
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In terms of Fourier transforms, we can write the above 
quantities as: 



n T /" 2r-V(cav). , s - 



Ej-j = lim 



dw 
2^' 



,2 T7(cav) 



-,^^cavJ/ / N - 

K (i-,r,w)e 



(C4) 



(C5) 



By the same arguments used earUer in the case 
A'^'^^^^ (r, r'), we see that Dij and Eij vanish, provided 

that ili'^^^^ and kI'^^^^ fall off like w~'^ or faster. By in- 
spection of Eqs. (jlOip and p02p . and recaUing that f(u') 
falls off like u>~^, we can easily see that this is the case 



for V 



(cav) 



Therefore we have: 



Consider now the quantity [7^''^^^'' whose expression is 
provided by Eq. (|100p . In order to estimate the fall- 
off rate of the various terms in the r.h.s. of Eq. (|100p . 
we need recall that Fresnel reflection coefficients of all 
real materials decay like w^"^ (see Appendix B). Since 
the quantities (1/^s — 1) and {l/Ap — 1) then decay 
like w^^, we see that all terms involving these quan- 
tities in the r.h.s. of Eq. (|100[) fall off at least like 
w^^ at large frequencies, and therefore they can be ne- 
glected. Consider now the remaining terms in the expres- 
sion for uI^'^^\y,y' ,w). Upon noticing that for large w, 



we get: 



k±ik±j + 0(w), 



(C6) 



o -rind 



(p) „ik(+'-r-ik<"'-r' , Jp) „ik("' ■r-ik(+' -r'+Sifc^ d 



k±ik±j + 0{w 



(C7) 



r 



Having reached this point, we take advantage of the fact 
that to order w^^ included, both Fresnel coefficients are 
independent of k± (see Eq. (|B2p V By virtue of this, the 
r..h.s. of Eq. (|C7p can be written as: 



(cav) 



^3 



^ E E Ki%lM+0{w-^) (C8) 



a—s,p m— 1,2 



where K^'^). . are material dependent constants, and 



m\ij 



Ii{w) — i 



2nk^ 



while 



h {w) = i 



2iTk 



gik<"'-r-ik<+>-r'-h2ifcj£i 



(C9) 



(CIO) 



also written as 

(fk 1 



and 



h{w) 



27r2 /fc2 _ fc2 



d^k 1 



(Cll) 



27r2 fc2 - fc2 



„j[ki-(ri-r^)+fc3(2<i-z-2')] 



(C12) 

From this we see that both Ii and I2 fall-off like w 2, 

V(cav) 



and therefore, in view of Eq. (jCSp . we find that t/,^ 
decays like w"^ or faster. Therefore Dij is zero, and then 
we obtain the desired result: 



clT^(r,r') = 0. 



(C13) 



It is a simple matter to check that the integrals h can be 
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